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Recently, the Fermilab heavy-quark action was extended to include dimension-six and -seven 
operators in order to reduce the discretization errors. In this talk, we present results of the first 
numerical simulations with this action (the OK action), where we study the masses of the quarko- 
nium and heavy-light systems. We calculate combinations of masses designed to test improve- 
ment and compare results obtained with the OK action to their counterparts obtained with the 
clover action. Our preliminary results show a clear improvement. 
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1. Introduction 

Simulating heavy quarks in lattice QCD is a challenging problem, because the quark mass niQ 
and the accessible ultraviolet cutoff a -1 are comparable. Special care is needed to handle dis- 
cretization errors [1]. In order to make accurate and reliable calculations of many Standard Model 
parameters involving heavy quarks, one needs not only computer power but also methodological 
improvements of the quark actions used. One line of attack is the Fermilab method [2], which 
stalls with the clover action [3] for Wilson fermions [4]. In the original work, interactions through 
dimension five were considered. More recently, we extended the Fermilab method to an action (the 
OK action) with dimension-six and -seven interactions [5, 6]. 

Using power counting as a guide, Ref. [5] estimated that the OK action should reduce dis- 
cretization effects for heavy quarks to ~ 1% on, say, the MILC asqtad ensembles [7]. In this paper, 
we present the first numerical results obtained with the OK action, to test whether the theoretical 
improvement is realized in practice. We compute combinations of rest masses and kinetic masses 
designed to test the improvement, without the need to tune the input parameters. 

In Sec. 2, we briefly discuss and present the OK action. Section 3 contains the details of the 
simulations and our preliminary results for an inconsistency / [8], and for the hyperfine splittings 
between rest and kinetic masses for the heavy-heavy and heavy-light systems. We discuss our 
results and future plans in section 4. 

2. OK Action and Tadpole Improvement 

In this section, we briefly describe the OK action, including the tadpole improvement [9] used 
in the simulation. In general, the Fermilab formalism calls for separate couplings for spatial and 
temporal interactions. Fermilab actions have a smooth transition to the m^a — > and m^a — > oo 
limits, but the short-distance coefficients depend on m^a, £ and the spatial Wilson parameter r s in 
a non-trivial way [2]. This action reduces the discretization errors to 0(a 2 A 2 b(mQa)), where b is 
a function that is bounded for all niQa. The OK action includes higher dimensional operators to 
further reduce the lattice spacing errors. 

Reference [5] starts by considering all the operators of dimension six and seven with two 
effective field theories in mind: heavy-quark effective theory (HQET) and the nonrelativistic QCD 
(NRQCD), appropriate to heavy-light and heavy-heavy systems, respectively. The interactions 
are classified in powers of A (HQET, X ~ A/rag, Aa) or v (NRQCD, relative internal velocity). 
Once all the independent operators are identified, the redundant ones are eliminated by means 
of field transformations, and the remaining couplings c; are determined via tree-level matching. 
Together with the one-loop matching of the dimension five chromomagnetic interaction, this action 
is expected to bring the discretization errors below the one-percent level [5]. 

For coding, and especially for tadpole improvement of the couplings, it is convenient to write 
the OK action in the hopping-parameter form: 

x x x 

- K t Y Vx[(r s C + 8c 4 ) - 7f(C - 2ci - l2c 2 ))TiXi/ x 

A. I 
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- *i£ %[(r£ + 8c 4 ) + 7l-(C " 2ci - 12c2)]r_i% 

jc,/ 

+ k> £ [2c 4 + (ci + 2c 2 )] 7] 2 y x - k, (c B C + 16c 5 ) £ y/ r /E • B M y x 

X X 

+ K t J^\j/ x {2c 4 -Yi(c 1 +2c2)]T* i y x -K t c E <;^W<x- E \atWx 

X X 

+ Ktc 2 £ ^•{7--r_ i ,r y -r_ j }v^+2K,c 5 £v^££{jX^ 

+ 2jqc 3 £ \j/ x {Y-D M ,iL • Bi a t}i^- + 2K- f c ££ £ v/i{7 4 £>4 lat , a • £i at }^v 

X X 

where Ti^y/^x) = U±^(x)\j/(x±afi) with U-^(x) = uj l {x — afl) and 

m a=- (l + 3r,C + 18c 4 ). (2.1) 

For further details and the matching conditions for the couplings c,, we refer the reader to Ref. [5]. 

Given the large one-loop corrections that can arise in lattice perturbation theory [9], before 
using the OK action in a numerical simulation we apply tadpole improvement to the matched 
couplings. To carry out the tadpole improvement, it is convenient to write T±^ = «o[7±/i/«o] = 
uqT±h where uq factors are absorbed into the couplings c,-. In this way, one finds the relations 
between bare and tadpole improved coefficients 





= U Q K t , 


(2.2) 


f,£ + 8c 4 


= r,C + 8c 4 


(2.3) 


£ — 2c\ — 12c2 


= £-2ci - \2c2 


(2.4) 


Q 




(2.5) 


c\ + 2c 2 


= u Q (a +2c 2 ) 


(2.6) 


C2 


= U C2, 


(2.7) 


c B £ + 16c 5 


= ul(c B £ + 16c 5 ) 


(2.8) 


CE 


3 

= UqCE, 


(2.9) 


C3 


= UqC3, 


(2.10) 


cee 


4 

— UqCee-, 


(2.11) 



where the last two follow because every term in the anticommutators has five links, while one 
power of Mo is absorbed, as usual, into k t . The matching conditions of Ref. [5] are then used, 
substituting 

moa = — - (1 + 3r,C + 18c 4 ) (2.12) 

for m$a. Another condition for c$ is needed, but it is not simple to express. In the expansion of B,i at 
in terms of Ts, one finds that vM^i-Bflat) (Tj + T-j)}Wx, J 7^ i. has terms with both 3 and 5 links. 
The 3-link terms arise from T±^T T ^ = 1. In the c 3 and cee interactions, each 3-link term appears 
twice, but with opposite sign, while here they have the same sign. Coding this operator with the 
correct uq factors is currently underway with USQCD software [10]. 



3 



Numerical Tests of the Improved Fermilab Action 



M. B. Oktay 



Table 1: Approximate values of the kinetic masses for the heavy-light and heavy-heavy systems obtained 
with the OK action. Thus, 0.042 is in the charm region and 0.036 is not yet in the bottom region. 





Mf> (Heavy-Light) [MeV] 


M| s (Heavy-Heavy) [MeV] 


0.036 


4418 


7332 


0.038 


3500 


5778 


0.040 


2680 


4227 


0.042 


1876 


2736 



3. Simulations and Tests 

We performed simulations on a "medium coarse' ' (a w 0.15 fm) 16 3 x 48 lattice with 2+1 
flavors of sea quarks, (ami,am s ) = (0.029,0.0484). From Ref. [11], we had data available with 
the clover action, and here we used the OK action with similar statistics, 500 configurations 
with 4 time sources per configuration. For the results presented here, we tadpole-improved the 
\ff x {iLjBj\ at , (Tj + T_j)}\jf x , j 7^ i, interaction with C5 = Mq c 5> pending completion of the code with 
proper tadpole improvement. We also choose f s = £ = £b = 1- This choice fixes the value of C4 
and C5 while the rest of the coefficients also depend on the choice of k t . 

Since the action is designed to improve 0(p 4 ) terms, we need to find observables to test these 
improvements. One such quantity is a combination of masses introduced in Ref. [8] and later 
discussed in Ref. [12]. Writing the rest and kinetic meson masses M^Q q and M 2 Q q as 

M \Qq = ™lQ+™lq+ B lQq, (3-1) 
M 2Qq = m 2Q+ m 2q+B 2 Qq, (3.2) 

where the ms are quark masses and the Bs are binding energies, Ref. [8] introduced the "inconsis- 
tency combination" 

i = 28M Qq - (SM QQ + 8M qq ) = 28B Qq - (8B Qq + 8B qq ) 
2M 2Qq 2M 2Qq 

where 8M = M2 — M\, 8B = B2 — B\ . The rightmost expression follows from the definitions. 

The binding energies B2 stem from the p 4 terms in the action [12]. By design, the the OK ac- 
tion improves these terms, compared with the clover action. Ideally, 8Bs and, hence, / should 
vanish, and we expect / to be smaller with an improved action. In order to compare the OK action 
with the clover action, we compute the rest and kinetic masses of heavy-light and heavy-heavy 
systems at four different hopping-parameter values, kt = 0.042,0.040,0.038,0.036. Table 1 lists 
the obtained pseudoscalar kinetic masses, to show the range of physical mass covered here. 

Figure 1 shows our results for the inconsistency /, together with results from an earlier study 
with the clover action [11]. For simplicity, and without serious loss in the strength of the test, we 
omit the light-light mass difference 8Mq q . As one can see, the clover action suffers from serious 
deviations from 1^0, especially for hopping parameters in the Z?-quark region. On the other 
hand, the OK action's inconsistency is statistically consistent with / = up to the largest masses 
considered. 
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Figure 1: Comparison of the inconsistency / between the OK action and the clover action. For the clover 
action K = 0.122 and K = 0.076 correspond to tuned rj 6 and r)/, masses. 



Another way to see the effects of improvement is to look at the hyperfine splittings. Let us 
define 

Ai = M^-M\ s , (3.4) 
A 2 =M^-M 2 PS , (3.5) 

where V and PS denote vector and pseudoscalar mesons. The rest-mass splitting Ai is accurate at 
the tree level, thanks to the clover term, but the kinetic-mass splitting A 2 has contributions from 
higher-dimension corrections such as {y-D,/E-B}, which are improved (unimproved) with the 
OK (clover) action. 

In Fig. 2 we plot aA 2 vs. aAi for quarkonium at the same values of k t as above. Ideally, the 
data would land on the line c?A 2 = aA] . We see that the OK data fare much better than the clover 
data. In Fig. 3 we plot aA 2 vs. aAi for a heavy-light meson. In this case, the clover action fares 
well to begin with, and our statistical errors are not small enough to test whether the OK action is 
an improvement. 



4. Outlook 



Our preliminary analysis of the OK action is encouraging. It shows clear improvements in 
the inconsistency / and the kinetic-mass hyperfine splittings in the heavy-heavy system. For the 
hyperfine splitting in the heavy-light system, the improvement is not yet clear, because the clover 
action already works well, so a decisive test requires higher statistics. The code for the OK inverter 
is still under development. The next step is to finish coding of the C5 term with the correct uq 
factors, to be followed by a thorough optimization. Thereafter, we plan on using the OK action for 
charm and bottom physics. 
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Figure 2: Kinetic -mass hyperfine splitting vs. rest-mass hyperfine splitting for the heavy-heavy system. 




25 I U lAJ I I I I 1 LJ LJ ' I I LJ ' I I I I ' U LJ ' I I I I ' LJ LJ ' I I LJ ' I ■ ■ ■ ' 

' 0.02 0.03 0.O4 0.05 0.06 0.07 0.08 0.09 0.1 0.11 0.12 

a[Mj(V) - M,(PS)] 



Figure 3: Kinetic-mass hyperfine splitting vs. rest-mass hyperfine splitting for the heavy-light system. 
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